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Abstract — High penetration levels of distributed photovoltaic 
(PV) generation on an electrical distribution circuit may severely 
degrade power quality due to voltage sags and swells caused by 
rapidly varying PV generation during cloud transients coupled 
with the slow response of existing utility compensation and 
regulation equipment. Although not permitted under current 
standards for interconnection of distributed generation, fast- 
reacting, VAR-capable PV inverters may provide the necessary 
reactive power injection or consumption to maintain voltage 
regulation under difficult transient conditions. As side benefit, 
the control of reactive power injection at each PV inverter 
provides an opportunity and a new tool for distribution utilities 
to optimize the performance of distribution circuits, e.g. by 
minimizing thermal losses. We suggest a local control scheme 
that dispatches reactive power from each PV inverter based on 
local instantaneous measurements of the real and reactive com- 
ponents of the consumed power and the real power generated by 
the PVs. Using one adjustable parameter per circuit, we balance 
the requirements on power quality and desire to minimize 
thermal losses. Numerical analysis of two exemplary systems, 
with comparable total PV generation albeit a different spatial 
distribution, show how to adjust the optimization parameter 
depending on the goal. Overall, this local scheme shows excellent 
performance; it's capable of guaranteeing acceptable power 
quality and achieving significant saving in thermal losses in 
various situations even when the renewable generation in excess 
of the circuit own load, i.e. feeding power back to the higher- 
level system. Key Words: Distributed Generation, Feeder Line, 
Power Flow, Voltage Control, Photovoltaic Power Generation 

I. Introduction 

The integration of significant quantities of renewable gen- 
eration into the U.S. electrical grids presents many significant 
challenges |1|. Reaching this goal will likely require a 
portfolio of renewable resources and generating methods, 
e.g. transmission-scale wind, concentrating solar power, and 
photovoltaic (PV) representing large scale options, while at 
the small or distribution scale, PV is presently the only 
viable option. Each of these different generation technologies 
constitutes a challenge which varies according to the location 
of its interconnection. At the transmission scale, renewable 
generation projects are generally large enough to warrant 
a transmission interconnection study. During the study, the 
transmission operator may uncover problems caused by the 
renewable generation and require the renewable generation 



owner to install certain equipment to mitigate the problems. 
In this case, the cost of mitigation is borne by the generator 
creating the transmission problem. 

At the distribution scale, the size of an individual PV gen- 
erator is so small that the cost of an "interconnection study" 
would be prohibitive. When the number of PV generators on 
a distribution circuit is small, the impact is quite small and 
the present utility systems are unaffected. However, as the PV 
penetration on a distribution circuit grows, the net impact of 
many small PV generators can reach a level where the power 
quality is significantly affected, e.g. loss voltage regulation 
due rapid variations in PV generation caused by cloud 
transients coupled with the slow response of existing utility 
equipment. Utilities could choose to install fast-response 
equipment to rectify the problem (e.g. a D-STATCOM El), 
however, the cost then would be socialized over the entire 
rate base and not borne by the PV generators creating the 
problem and benefiting economically from interconnection 
to the distribution grid. 

An alternate solution is to place the burden of providing 
reactive power for voltage regulation on the individual, small- 
scale PV generators by using excess PV inverter capacity 
to generate or consume reactive power. However, purposeful 
injection of reactive power or attempting to regulate volt- 
age by a distributed generator is not permitted by current 
interconnection standards 0. Changes to these standards 
to allow for controlling reactive power appear eminent, but 
even if these changes are made, there still remains significant 
questions as to how much reactive power to dispatch from 
each inverter, when to dispatch it, and where and how the 
control signals are generated. 

In contrast to the transmission system, the number and 
diversity of distribution circuits is perhaps too large to model 
the individual placement of PV generators. Therefore, to 
assess the effectiveness of different algorithms for dispatch- 
ing reactive power from PV inverters, we take a statistical 
sampling approach. We fix certain average 'macroscopic' 
parameters of the distribution circuit: spacing between loads, 
real and reactive power drawn per load, and PV generation. 
However, when creating a realization of the circuit, the 



individual characteristics of links and values of loads are 
selected at random from assumed distributions. For a given 
set of macroscopic parameters, we create sample realiza- 
tions of load and generation profiles and analyze the circuit 
response under the reactive power dispatch algorithm. The 
macroscopic parameters are varied to assess many different 
scenarios, e.g. when there are many small PV generators on 
the circuit but their output is low compared to the circuit 
loading and when there are only a few PV panels on the 
feeder line but they inject a substantial amount of power in 
the circuit leading a net power export and a rise of local 
voltage. 

Analysis such as described above can assess the per- 
formance of different algorithms over a wide variety of 
distribution circuits and can provide valuable guidance to 
regulatory bodies on questions such as how much excess PV 
inverter capacity is required to ensure voltage regulation on 
wide array of distribution circuits. 

In our recent work, [4 ] we have focused on the comparison 
of centralized and decentralized (local) approaches to the 
control of reactive power. We have shown that, for realistic 
feeder lines, a simple local control technique can achieve 
almost 80% of savings in losses when compared to a central- 
ized control based on solving the full optimization problem. 
In this work, we extend this approach by discussing the effect 
of a class of local control techniques on both the losses and 
power quality in the system. Although our particular reactive 
power dispatch algorithm may not work in all cases, the main 
purpose of this manuscript is to develop a framework for 
analyzing the effectiveness of any algorithm. 

There are other approaches to optimizing the dispatch of 
reactive power for the purpose of voltage regulation and loss 
minimization that could be adapted to the present problem 
including work by Baran and Wu0, 0, (7) and Baldick 
and Wu (8) and also in (9), ifTUl . However, these works are 
somewhat specialized to optimal placement and/or control a 
few large sources of reactive power where the problem at 
hand includes many small sources. 

The remainder of this manuscript is organized as follows. 
Section [II] describes the capability of an inverter to inject or 
consume reactive power. Section [TIT] describes the power flow 
solution method we adopt and introduces our reactive power 



dispatch algorithm. Section [IV] provides a description of the 
prototypical distribution circuit where we test our algorithm. 
Section [V] presents the results of our simulations and shows 
under what conditions our reactive power dispatch algorithm 



performs well, and Section VI gives our conclusions and 
suggestions for future work. 

II. Inverter as a limited regulator of local 

REACTIVE POWER FLOW. 

Under the interconnection standard for distributed gener- 
ation [3], PV inverters may not inject or consume reactive 
power or attempt to regulate voltage in any way, i.e. they must 
operate at unity power factor while in grid-tied mode. To 
overcome voltage regulation problems anticipated on distri- 
bution circuits with high penetrations of PV generation, these 
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Fig. 1. When s is larger than p( g \ the inverter can supply or consume 
reactive power q^ 9 \ The inverter can dispatch q^) quickly (on the cycle- 
to-cycle time scale) providing a mechanism for rapid voltage regulation. As 
the output of the PV panel array p( g ) approaches s, the range of available 
g^) decreases to zero. 



regulations are expected to be changed to allow the injection 
of reactive power. Before we discuss how to dispatch reactive 
power and analyze the effectiveness of the dispatch, we first 
discuss the limitations on a PV inverter's reactive power 
capability. 

We adopt a model of PV inverters previously described in 
HI . ifTTIl . If the apparent power capability s of an inverter 
exceed the instantaneous real power generated p^ by the 
connected PV panels, the range of allowable reactive power 
generation is given by \q^\ < \/ s 2 — {p^)) 2 = q max . This 
relationship is also described by the phasor diagram in Fig.[T] 
On a clear day with the sun angle aligned with the PV array, 
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and the range of available is at a minimum. 
One focus of this manuscript is to provide a framework for 
statistically analyzing distribution circuits so set a minimum 
s relative to pmax- However, the analysis is complicated 
because during cloud events or when the sun angle is not 



perfectly aligned with the PV panels, p^ 9) < pmax implying 
that the range of reactive power capability varies though the 
year, the day, and with the weather. In order to model this 
effect we perform our simulations for different values of p^ 
while keeping the the absolute capability s fixed. 

III. Optimization of Losses and Voltage Control 

To solve for complex power flows in a radial circuit, we 
adopt the DistFlow AC power flow equations of Q, (6), 
Q, which for the radial circuit illustrated in Fig. ^ can be 
written Vj = 1, • • • , n: 
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Fig. 2. Diagram and notations for the radial network. Pj and Qj represent 
real and reactive power flowing down the circuit from node j, where Prj and 
Qo represent the power flow from the sub-station, pj and qj correspond to 
the flow of power out of the network at the node j, where the respective 
positive [negative] contributions, p^y* and \Pj and q^] represent 
consumption [generation] of power at the node. The node-local control 
parameter q^ can be positive or negative but is bounded in absolute value 
as described in Eq. [4] The apparent power capability of the inverter Sj is 



(<?) 



preset to a value comparable to but larger than maxpV 



where Pj + iQj is the complex power flowing away from 
node j toward node j + 1, Vj is the voltage at node j, Tj +%Xj 
is the complex impedance of the link between node j and 
j + 1, and pj +iqj is the complex power extracted at the node 
j. Both pj and qj are composed of local consumption minus 
local generation due to the PV inverter, i.e. pj = p^ — p^ 
and qj = q^ — q^\ Of the four contributions to pj + iqj, 

Pj C ^ , Pj 9 ^ and, q^ are uncontrolled (i.e. driven by consumer 
load or instantaneous PV generation), while the reactive 
power generated by the PV inverter, q^\ can be adjusted. 

As described in Section [n| q^ is limited by the reactive 
capability of the inverter: 



Vj = 1, 



< 



(4) 



Note that reactive power generation is possible only at the 
nodes with PV generation. For the other nodes, we take sj — 
0. For a broad class of feeder lines such as the one considered 
in this work, the quadratic terms in Eqs. ( 1|2|3 } are relatively 
small as are the deviations of voltage along the line. In this 
case the power flow equations can be approximated by the 
linear LinDistFlow equations 0, 0, Q: 
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where we have exploited the approximation V fe 2 
2Vb (V k - V ). 

Within the framework of linearized model, the rate of 
energy dissipation (losses) in the distribution circuit is given 
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Minimizing or at least keeping the losses acceptably low is a 
natural goal for optimization and control. However, voltage 
variations along the circuit must stay within strict regulation 
bounds. These bounds are respected if the maximum devia- 
tion of the per unit voltage SV obeys 
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where e w 0.05 in normal operation. The goal of the voltage 
stability optimization is to keep SV within normal bounds. 

In this work, we discuss local control techniques that try 
to achieve objectives of minimizing both the thermal losses 
^ and the maximal voltage deviation ([9]). Our local scheme 
should be viewed as an approximate/heuristic solution of the 
following multi-objective optimization problem 



min [£, SV] T , 
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s.t.Eqs. ( |4|5|6|7| ). 



The locality of the control scheme corresponds to the re- 
quirement that the q^ depend only on local information, 
which we restrict to p^\ p^\ qjf* . Another local variable 
that could be used is Vk, however, using Vk could easily 
lead to inequities. For instance, compared to a PV generator 
near the beginning of a distribution circuit, one at the end 
of the circuit may perpetually see a relatively lower voltage 
and be requested to inject reactive power leading to faster 
degradation of the inverter. Instead, our scheme is based on 
local real and reactive power consumption and generation 
and should be more equitable along the circuit. We also 
assume that the control scheme is homogeneous over the line: 
all inverters are programmed in the same way, and explicit 
dependence on the bus number k enters through the inverter 
capability s^. Formally, we will study local control schemes 
with 
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and consistent with constraint ([4}. It is useful to introduce 
the following "helper" function, Constr^, meant to enforce 
the constraint ([4]): 
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The local control scheme from (H was designed based on the 
idea that losses are minimized when the reactive flows Qk 
are zero. The qjf^ were chosen to minimize the net reactive 
power consumption q^ 
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at each node: 
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This scheme was shown to be very effective in reducing the 
losses. In this work, we extend this idea and propose a new 
class of the schemes that attempt to achieve simultaneously 
both optimization objectives in Eq. ( [T0| ). 



Eq. ([7} suggests that, to reduce variations in we should 
minimize the absolute value of the combined power flow 
VkPu + XkQk- Note that for many circuits, the ratio of 
Tk/xk = ol is nearly constant for all k and depends only on 
the configuration and size of the conductors used. Thus the 
absolute value of VkPk+XkQk will be exactly zero if for ev- 
ery line we will ensure that p^ —p^ +a — q^^j = 0. 
Summarizing, this suggests the following control function 
F^ V \ aimed at minimizing voltage variations and ignoring 
losses: 



In essence, is attempting to not only supply the reactive 
power needs of the local load at node k, but also to supply 
the reactive power consumption of the adjacent links of 
the distribution circuit. A compromise between the two 
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objectives in ( 12 ) can be achieved via the following nonlinear 
combination 



F k = Constr/e 



KF^ + {l-K)F } 
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where K is a single parameter controlling the trade off 



between the two objectives in Eq. ( 12 ). At K — 1 we recover 
the fl3] ) scheme, whereas at K = the scheme reduces to 

IV. Description of the prototypical rural 

DISTRIBUTION CIRCUIT 

We consider a sparsely-loaded rural distribution circuit 
model with 250 nodes based on one of the 24 prototypical 
distribution circuits described in fT2l . The nominal phase-to- 
neutral voltage Vq is 1.2kV and the line impedance is set 
to (0.33 + 0.38i)Q/km, constant along the distribution line. 
The distance between two neighboring nodes is uniformly 
distributed between 0.2 and 0.3 kilometers. The real powers 

(c) 

Pj consumed by the loads are uniformly distributed between 



and p 



) where p( c > max ) is set to either to lkW or 



is 



2.bkW (for the different cases we consider), and q 
randomly selected from uniform distribution between 0.2p 

(c) 

and 0.3p!- ' . In this work, we study two levels of PV pene- 
tration; either 20% or 50% of the nodes have PV generation 
installed, but the PV-enabled nodes are selected randomly. We 
assume that the generated power p^ is the same for all PV- 
enabled nodes being equal to either lkW or2kW. A uniform 
level of PV power generation assumes identical installations 
at each PV-enabled node (the same pmax installed in the same 
way) and uniform solar irradiance. The inverter capacity is 
set to s = 2.2 k W at PV-enabled nodes. 

The algorithm we use for the simulation consists of several 
steps: first we generate the random values of r& and Xk that 
remain constant for all forthcoming simulations. For every 
case, we generate a random sample of loads and generation 
and solve the LinDistFlow equations ([5][7]) to find both the 
voltage levels along the line and the total losses. We do not 
perform any averaging over many realizations of the load 
distribution because the results do not differ significantly 



TABLE I 

Comparison of parameters for different model feeder lines 



Case # 


PR 


p(c,mean) 


p {9) 


SV 


Co 


1 


20 % 


1.25 kW 


1.0 kW 


0.059 


7.84 kW 


2 


20 % 


0.5 kW 


2.0 kW 


0.014 


0.33 kW 


3 


50 % 


1.25 kW 


1.0 kW 


0.048 


4.66 kW 


4 


50 % 


0.5 kW 


2.0 kW 


0.014 


1.89 kW 



Case 1 
Case 2 
Case 3 
Case 4 



Fig. 3. 



SV 



max/j 



Vfc-Vo I 



K 



vs K for the local control according to 



fT5) and four cases described in Table [I] 



from one sample to another due to the self- averaging effect 
of a large number of node on the line. 

V. Simulations: Results and Discussions 

Four regimes of operation are considered corresponding to 
two different levels of PV penetration and two operational 
extremes (a) sunny day, everybody is at work, and (b) 
overcast, everybody home. These are summarized in Table 
[I] In cases 1 and 3 the amount of power generated is 
low compared to the total load of the system, however the 
capacity of the inverters is high: s/p^ « 2.2. In the other 
two cases the situation is opposite; the total amount of power 
generated is almost equal to the total consumption in case 2, 
and there is over generation of power in case 4 leading to 
reversal of the power flow (the circuit is feeding power into 
the higher level grid). The inverter capacity in cases 2 and 4 
is relatively small: s/p^ « 1.1. These 4 cases allow us to 
assess the importance of several characteristics of the system: 
generation and load levels, PV penetration rate, and inverter 
capacity. 



The reactive power is controlled according to Eq. ( [15] ) with 
the parameter K set to the same level on all the generating 
nodes. K is varied in the range — 5 < K < 10, and the total 
losses C and maximal voltage deviation SV are computed. 
For all cases, the performance of the control technique is 
compared to the base case of zero reactive power generation: 
q^ =0. The corresponding values of Co and SVo are shown 
in Table U 

The dependence 5V versus K is shown in Fig. [3] The 
performance of the voltage control is very good in cases 2 




Fig. 4. Dissipated power with control, measured in the units of the 
respective dissipated power without control, vs K for the local control 
according to (l"5) and four regimes described in Table [I] The lowest relative 
losses achievable are 0.93, 0.67, 0.86 and 0.94 for cases 1 to 4 respectively. 



and 4, when there is a small amount of power consumption. 
Significant reduction in voltage deviation in the case 4 is 
possible because of almost complete compensation of the 
reversed flow of real power by increase in the consumption 
of reactive power. In cases 1 and 3, correspondent to high 
power consumption, the reduction in voltage deviation is less 
significant, however it is still possible to achieve the reduction 
of ^0.01 in both cases. 

Changes in losses with variation in the parameter, K, 
shown in Fig. [4] are similar to these of the power quality 
discussed above. The top savings are achieved in the cases 

2 and 3. Thus, at the optimal point of the case 2 the con- 
sumption is almost completely balanced by the generation, 
and then the magnitudes of total power flows Pk-,Qk are 
greatly reduced. Note, that although the relative savings in 
losses are very impressive in cases 2 and 4, the absolute 
values are still smaller compared to these in the cases 1 and 

3 because of the lower values of Co in the latter. The high 
absolute savings in the cases 1 and 3 should be attributed to 
the high value of the relative inverter capacity: s/p^ =2.2, 
that allows high injections of reactive power in the system. 
Despite, the high penetration level the savings in losses in 
the case 4 are less pronounced in comparison to the case 2 
because of significant losses associated with the reverse flow 
of real power. Interestingly, and in a contrast to other cases, 
optimization of losses in the case 4 required fine tuning of 
the parameter K. 

The overall performance of the control scheme [T5]can be 
most easily assessed on the combined plot on Figure|5] where 
the efficiency in both SV and C is shown simultaneously. 
It is apparent from this plot that in the cases 1 and 3 the 
optimum in both losses and voltage deviations is achieved 
simultaneously. From the viewpoint of the multi-objective 
optimization this corresponds to Pareto set consisting of a 
single point. The situation is different in the cases 2 and 4 
where the Pareto set is represented on a range in K. (This 
set corresponds to the negative slope parts of the curves in 



Fig. 5. The plot combines results from Fig.|4]and Fig.[3]in one, thus showing 
dependence of the relative power loss vs power quality characteristic, when 
scanned over allowed range of the parameter K in control scheme^] In 
the cases 1 and 3, the left most point of the curve corresponds to the 
best configuration, achieving both the lowest 5V and the smallest losses. 
For cases 3 and 4, the whole part of the curve from the left most point 
(corresponding to the best power quality) to the point where global minimum 
is achieved (losses are minimal then) represents a feasible region where one 
can trade-off between savings in losses and the voltage drop by adjusting 
K. 



Fig. [5]) Comparing different cases, the maximal flexibility 
is achieved in the case 4 where the high penetration of 
renewables results in the overgeneration of power. 

VI. Conclusions and Path Forward 



To summarize, this study suggested 



• Simple and local control scheme (15) for the reactive 
power generation at circuit-distributed PV inverters. The 
control adjusts the reactive power according to the local 
values of (real and reactive) consumption and real PV 
generation. The scheme contains one global adjustable 
parameter balancing between complementary (local) de- 
sires to minimize the power flow and to maintain good 
power quality. 

• We experimented with the scheme showing that the 
local scheme is capable of achieving very significant 
(and probably sufficiently close to globally optimal) 
simultaneous improvement in global power quality and 
global reduction of losses over the entire circuit. Dif- 
ferent regimes tested correspond to sunny/overcast and 
large/small load cases. 

• We concluded that already 10% of excess inverter ca- 
pacity s is enough to allow significant reductions in both 
the losses and voltage deviations. 

• We observe that in the case of the over-generation, and 
thus reversal of the power flow direction, the multi- 
objective optimization becomes more sensitive to the 
parameter change, thus requiring a more accurate tuning. 

The results reported in this study will require further 
analysis and exploration along the following lines: 

• Study of other local control schemes aimed at even 
better improvements. 



• We conjectured that fluctuations from an instance to 
instance (e.g. in configurations of load) drawn from 
a reasonable distribution representing a typical regime 
(say a morning of the overcast day) will not be sig- 
nificant. This conjecture needs to be tested in a more 
accurate statistical study. The study should test statistical 
variations in the control scheme performance. 

• It may be beneficial for improving the power quality 
characteristics to include information about local volt- 
age (if such information is available) into the control 
scheme. 

• We also envision generalizing the scheme to a global 
dynamical control, with the parameter K varying in time 
according to a broadcasted central command. Moreover, 
a distributed control, with the coefficient K also varying 
spatially, may also be considered as beneficial in special 
(geographically varying) situations, where we should 
also consider inhomogeneity and significant variation in 
the instantaneous insulation pattern. 
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